3.7.5 \(\int x^4 (d+e x^2)^2 (a+b \arcsin (c x)) \, dx\) [605]

3.7.5.1 Optimal result
3.7.5.2 Mathematica [A] (verified)
3.7.5.3 Rubi [A] (verified)
3.7.5.4 Maple [A] (verified)
3.7.5.5 Fricas [A] (verification not implemented)
3.7.5.6 Sympy [A] (verification not implemented)
3.7.5.7 Maxima [A] (verification not implemented)
3.7.5.8 Giac [B] (verification not implemented)
3.7.5.9 Mupad [F(-1)]

3.7.5.1 Optimal result

Integrand size = 21, antiderivative size = 241 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\frac {b \left (63 c^4 d^2+90 c^2 d e+35 e^2\right ) \sqrt {1-c^2 x^2}}{315 c^9}-\frac {2 b \left (63 c^4 d^2+135 c^2 d e+70 e^2\right ) \left (1-c^2 x^2\right )^{3/2}}{945 c^9}+\frac {b \left (21 c^4 d^2+90 c^2 d e+70 e^2\right ) \left (1-c^2 x^2\right )^{5/2}}{525 c^9}-\frac {2 b e \left (9 c^2 d+14 e\right ) \left (1-c^2 x^2\right )^{7/2}}{441 c^9}+\frac {b e^2 \left (1-c^2 x^2\right )^{9/2}}{81 c^9}+\frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x)) \]

output
-2/945*b*(63*c^4*d^2+135*c^2*d*e+70*e^2)*(-c^2*x^2+1)^(3/2)/c^9+1/525*b*(2 
1*c^4*d^2+90*c^2*d*e+70*e^2)*(-c^2*x^2+1)^(5/2)/c^9-2/441*b*e*(9*c^2*d+14* 
e)*(-c^2*x^2+1)^(7/2)/c^9+1/81*b*e^2*(-c^2*x^2+1)^(9/2)/c^9+1/5*d^2*x^5*(a 
+b*arcsin(c*x))+2/7*d*e*x^7*(a+b*arcsin(c*x))+1/9*e^2*x^9*(a+b*arcsin(c*x) 
)+1/315*b*(63*c^4*d^2+90*c^2*d*e+35*e^2)*(-c^2*x^2+1)^(1/2)/c^9
 
3.7.5.2 Mathematica [A] (verified)

Time = 0.13 (sec) , antiderivative size = 187, normalized size of antiderivative = 0.78 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\frac {315 a x^5 \left (63 d^2+90 d e x^2+35 e^2 x^4\right )+\frac {b \sqrt {1-c^2 x^2} \left (4480 e^2+160 c^2 e \left (81 d+14 e x^2\right )+24 c^4 \left (441 d^2+270 d e x^2+70 e^2 x^4\right )+4 c^6 \left (1323 d^2 x^2+1215 d e x^4+350 e^2 x^6\right )+c^8 \left (3969 d^2 x^4+4050 d e x^6+1225 e^2 x^8\right )\right )}{c^9}+315 b x^5 \left (63 d^2+90 d e x^2+35 e^2 x^4\right ) \arcsin (c x)}{99225} \]

input
Integrate[x^4*(d + e*x^2)^2*(a + b*ArcSin[c*x]),x]
 
output
(315*a*x^5*(63*d^2 + 90*d*e*x^2 + 35*e^2*x^4) + (b*Sqrt[1 - c^2*x^2]*(4480 
*e^2 + 160*c^2*e*(81*d + 14*e*x^2) + 24*c^4*(441*d^2 + 270*d*e*x^2 + 70*e^ 
2*x^4) + 4*c^6*(1323*d^2*x^2 + 1215*d*e*x^4 + 350*e^2*x^6) + c^8*(3969*d^2 
*x^4 + 4050*d*e*x^6 + 1225*e^2*x^8)))/c^9 + 315*b*x^5*(63*d^2 + 90*d*e*x^2 
 + 35*e^2*x^4)*ArcSin[c*x])/99225
 
3.7.5.3 Rubi [A] (verified)

Time = 0.47 (sec) , antiderivative size = 241, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.238, Rules used = {5230, 27, 1578, 1195, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx\)

\(\Big \downarrow \) 5230

\(\displaystyle -b c \int \frac {x^5 \left (35 e^2 x^4+90 d e x^2+63 d^2\right )}{315 \sqrt {1-c^2 x^2}}dx+\frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 27

\(\displaystyle -\frac {1}{315} b c \int \frac {x^5 \left (35 e^2 x^4+90 d e x^2+63 d^2\right )}{\sqrt {1-c^2 x^2}}dx+\frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 1578

\(\displaystyle -\frac {1}{630} b c \int \frac {x^4 \left (35 e^2 x^4+90 d e x^2+63 d^2\right )}{\sqrt {1-c^2 x^2}}dx^2+\frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 1195

\(\displaystyle -\frac {1}{630} b c \int \left (\frac {35 e^2 \left (1-c^2 x^2\right )^{7/2}}{c^8}-\frac {10 e \left (9 d c^2+14 e\right ) \left (1-c^2 x^2\right )^{5/2}}{c^8}+\frac {3 \left (21 d^2 c^4+90 d e c^2+70 e^2\right ) \left (1-c^2 x^2\right )^{3/2}}{c^8}-\frac {2 \left (63 d^2 c^4+135 d e c^2+70 e^2\right ) \sqrt {1-c^2 x^2}}{c^8}+\frac {63 d^2 c^4+90 d e c^2+35 e^2}{c^8 \sqrt {1-c^2 x^2}}\right )dx^2+\frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{5} d^2 x^5 (a+b \arcsin (c x))+\frac {2}{7} d e x^7 (a+b \arcsin (c x))+\frac {1}{9} e^2 x^9 (a+b \arcsin (c x))-\frac {1}{630} b c \left (\frac {20 e \left (1-c^2 x^2\right )^{7/2} \left (9 c^2 d+14 e\right )}{7 c^{10}}-\frac {70 e^2 \left (1-c^2 x^2\right )^{9/2}}{9 c^{10}}-\frac {6 \left (1-c^2 x^2\right )^{5/2} \left (21 c^4 d^2+90 c^2 d e+70 e^2\right )}{5 c^{10}}+\frac {4 \left (1-c^2 x^2\right )^{3/2} \left (63 c^4 d^2+135 c^2 d e+70 e^2\right )}{3 c^{10}}-\frac {2 \sqrt {1-c^2 x^2} \left (63 c^4 d^2+90 c^2 d e+35 e^2\right )}{c^{10}}\right )\)

input
Int[x^4*(d + e*x^2)^2*(a + b*ArcSin[c*x]),x]
 
output
-1/630*(b*c*((-2*(63*c^4*d^2 + 90*c^2*d*e + 35*e^2)*Sqrt[1 - c^2*x^2])/c^1 
0 + (4*(63*c^4*d^2 + 135*c^2*d*e + 70*e^2)*(1 - c^2*x^2)^(3/2))/(3*c^10) - 
 (6*(21*c^4*d^2 + 90*c^2*d*e + 70*e^2)*(1 - c^2*x^2)^(5/2))/(5*c^10) + (20 
*e*(9*c^2*d + 14*e)*(1 - c^2*x^2)^(7/2))/(7*c^10) - (70*e^2*(1 - c^2*x^2)^ 
(9/2))/(9*c^10))) + (d^2*x^5*(a + b*ArcSin[c*x]))/5 + (2*d*e*x^7*(a + b*Ar 
cSin[c*x]))/7 + (e^2*x^9*(a + b*ArcSin[c*x]))/9
 

3.7.5.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 1195
Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))^(n_.)*((a_.) + (b_.)*(x 
_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d + e*x)^m*(f + 
 g*x)^n*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, n}, x 
] && IGtQ[p, 0]
 

rule 1578
Int[(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_ 
)^4)^(p_.), x_Symbol] :> Simp[1/2   Subst[Int[x^((m - 1)/2)*(d + e*x)^q*(a 
+ b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x] && Int 
egerQ[(m - 1)/2]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 5230
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_ 
)^2)^(p_.), x_Symbol] :> With[{u = IntHide[(f*x)^m*(d + e*x^2)^p, x]}, Simp 
[(a + b*ArcSin[c*x])   u, x] - Simp[b*c   Int[SimplifyIntegrand[u/Sqrt[1 - 
c^2*x^2], x], x], x]] /; FreeQ[{a, b, c, d, e, f, m}, x] && NeQ[c^2*d + e, 
0] && IntegerQ[p] && (GtQ[p, 0] || (IGtQ[(m - 1)/2, 0] && LeQ[m + p, 0]))
 
3.7.5.4 Maple [A] (verified)

Time = 0.47 (sec) , antiderivative size = 329, normalized size of antiderivative = 1.37

method result size
parts \(a \left (\frac {1}{9} e^{2} x^{9}+\frac {2}{7} d e \,x^{7}+\frac {1}{5} d^{2} x^{5}\right )+\frac {b \left (\frac {c^{5} \arcsin \left (c x \right ) e^{2} x^{9}}{9}+\frac {2 c^{5} \arcsin \left (c x \right ) e d \,x^{7}}{7}+\frac {\arcsin \left (c x \right ) d^{2} c^{5} x^{5}}{5}-\frac {35 e^{2} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )+63 d^{2} c^{4} \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )+90 d \,c^{2} e \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )}{315 c^{4}}\right )}{c^{5}}\) \(329\)
derivativedivides \(\frac {\frac {a \left (\frac {1}{5} d^{2} c^{9} x^{5}+\frac {2}{7} d \,c^{9} e \,x^{7}+\frac {1}{9} e^{2} c^{9} x^{9}\right )}{c^{4}}+\frac {b \left (\frac {\arcsin \left (c x \right ) d^{2} c^{9} x^{5}}{5}+\frac {2 \arcsin \left (c x \right ) d \,c^{9} e \,x^{7}}{7}+\frac {\arcsin \left (c x \right ) e^{2} c^{9} x^{9}}{9}-\frac {e^{2} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )}{9}-\frac {d^{2} c^{4} \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )}{5}-\frac {2 d \,c^{2} e \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )}{7}\right )}{c^{4}}}{c^{5}}\) \(339\)
default \(\frac {\frac {a \left (\frac {1}{5} d^{2} c^{9} x^{5}+\frac {2}{7} d \,c^{9} e \,x^{7}+\frac {1}{9} e^{2} c^{9} x^{9}\right )}{c^{4}}+\frac {b \left (\frac {\arcsin \left (c x \right ) d^{2} c^{9} x^{5}}{5}+\frac {2 \arcsin \left (c x \right ) d \,c^{9} e \,x^{7}}{7}+\frac {\arcsin \left (c x \right ) e^{2} c^{9} x^{9}}{9}-\frac {e^{2} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )}{9}-\frac {d^{2} c^{4} \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )}{5}-\frac {2 d \,c^{2} e \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )}{7}\right )}{c^{4}}}{c^{5}}\) \(339\)

input
int(x^4*(e*x^2+d)^2*(a+b*arcsin(c*x)),x,method=_RETURNVERBOSE)
 
output
a*(1/9*e^2*x^9+2/7*d*e*x^7+1/5*d^2*x^5)+b/c^5*(1/9*c^5*arcsin(c*x)*e^2*x^9 
+2/7*c^5*arcsin(c*x)*e*d*x^7+1/5*arcsin(c*x)*d^2*c^5*x^5-1/315/c^4*(35*e^2 
*(-1/9*c^8*x^8*(-c^2*x^2+1)^(1/2)-8/63*c^6*x^6*(-c^2*x^2+1)^(1/2)-16/105*c 
^4*x^4*(-c^2*x^2+1)^(1/2)-64/315*c^2*x^2*(-c^2*x^2+1)^(1/2)-128/315*(-c^2* 
x^2+1)^(1/2))+63*d^2*c^4*(-1/5*c^4*x^4*(-c^2*x^2+1)^(1/2)-4/15*c^2*x^2*(-c 
^2*x^2+1)^(1/2)-8/15*(-c^2*x^2+1)^(1/2))+90*d*c^2*e*(-1/7*c^6*x^6*(-c^2*x^ 
2+1)^(1/2)-6/35*c^4*x^4*(-c^2*x^2+1)^(1/2)-8/35*c^2*x^2*(-c^2*x^2+1)^(1/2) 
-16/35*(-c^2*x^2+1)^(1/2))))
 
3.7.5.5 Fricas [A] (verification not implemented)

Time = 0.26 (sec) , antiderivative size = 219, normalized size of antiderivative = 0.91 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\frac {11025 \, a c^{9} e^{2} x^{9} + 28350 \, a c^{9} d e x^{7} + 19845 \, a c^{9} d^{2} x^{5} + 315 \, {\left (35 \, b c^{9} e^{2} x^{9} + 90 \, b c^{9} d e x^{7} + 63 \, b c^{9} d^{2} x^{5}\right )} \arcsin \left (c x\right ) + {\left (1225 \, b c^{8} e^{2} x^{8} + 10584 \, b c^{4} d^{2} + 50 \, {\left (81 \, b c^{8} d e + 28 \, b c^{6} e^{2}\right )} x^{6} + 12960 \, b c^{2} d e + 3 \, {\left (1323 \, b c^{8} d^{2} + 1620 \, b c^{6} d e + 560 \, b c^{4} e^{2}\right )} x^{4} + 4480 \, b e^{2} + 4 \, {\left (1323 \, b c^{6} d^{2} + 1620 \, b c^{4} d e + 560 \, b c^{2} e^{2}\right )} x^{2}\right )} \sqrt {-c^{2} x^{2} + 1}}{99225 \, c^{9}} \]

input
integrate(x^4*(e*x^2+d)^2*(a+b*arcsin(c*x)),x, algorithm="fricas")
 
output
1/99225*(11025*a*c^9*e^2*x^9 + 28350*a*c^9*d*e*x^7 + 19845*a*c^9*d^2*x^5 + 
 315*(35*b*c^9*e^2*x^9 + 90*b*c^9*d*e*x^7 + 63*b*c^9*d^2*x^5)*arcsin(c*x) 
+ (1225*b*c^8*e^2*x^8 + 10584*b*c^4*d^2 + 50*(81*b*c^8*d*e + 28*b*c^6*e^2) 
*x^6 + 12960*b*c^2*d*e + 3*(1323*b*c^8*d^2 + 1620*b*c^6*d*e + 560*b*c^4*e^ 
2)*x^4 + 4480*b*e^2 + 4*(1323*b*c^6*d^2 + 1620*b*c^4*d*e + 560*b*c^2*e^2)* 
x^2)*sqrt(-c^2*x^2 + 1))/c^9
 
3.7.5.6 Sympy [A] (verification not implemented)

Time = 1.30 (sec) , antiderivative size = 415, normalized size of antiderivative = 1.72 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\begin {cases} \frac {a d^{2} x^{5}}{5} + \frac {2 a d e x^{7}}{7} + \frac {a e^{2} x^{9}}{9} + \frac {b d^{2} x^{5} \operatorname {asin}{\left (c x \right )}}{5} + \frac {2 b d e x^{7} \operatorname {asin}{\left (c x \right )}}{7} + \frac {b e^{2} x^{9} \operatorname {asin}{\left (c x \right )}}{9} + \frac {b d^{2} x^{4} \sqrt {- c^{2} x^{2} + 1}}{25 c} + \frac {2 b d e x^{6} \sqrt {- c^{2} x^{2} + 1}}{49 c} + \frac {b e^{2} x^{8} \sqrt {- c^{2} x^{2} + 1}}{81 c} + \frac {4 b d^{2} x^{2} \sqrt {- c^{2} x^{2} + 1}}{75 c^{3}} + \frac {12 b d e x^{4} \sqrt {- c^{2} x^{2} + 1}}{245 c^{3}} + \frac {8 b e^{2} x^{6} \sqrt {- c^{2} x^{2} + 1}}{567 c^{3}} + \frac {8 b d^{2} \sqrt {- c^{2} x^{2} + 1}}{75 c^{5}} + \frac {16 b d e x^{2} \sqrt {- c^{2} x^{2} + 1}}{245 c^{5}} + \frac {16 b e^{2} x^{4} \sqrt {- c^{2} x^{2} + 1}}{945 c^{5}} + \frac {32 b d e \sqrt {- c^{2} x^{2} + 1}}{245 c^{7}} + \frac {64 b e^{2} x^{2} \sqrt {- c^{2} x^{2} + 1}}{2835 c^{7}} + \frac {128 b e^{2} \sqrt {- c^{2} x^{2} + 1}}{2835 c^{9}} & \text {for}\: c \neq 0 \\a \left (\frac {d^{2} x^{5}}{5} + \frac {2 d e x^{7}}{7} + \frac {e^{2} x^{9}}{9}\right ) & \text {otherwise} \end {cases} \]

input
integrate(x**4*(e*x**2+d)**2*(a+b*asin(c*x)),x)
 
output
Piecewise((a*d**2*x**5/5 + 2*a*d*e*x**7/7 + a*e**2*x**9/9 + b*d**2*x**5*as 
in(c*x)/5 + 2*b*d*e*x**7*asin(c*x)/7 + b*e**2*x**9*asin(c*x)/9 + b*d**2*x* 
*4*sqrt(-c**2*x**2 + 1)/(25*c) + 2*b*d*e*x**6*sqrt(-c**2*x**2 + 1)/(49*c) 
+ b*e**2*x**8*sqrt(-c**2*x**2 + 1)/(81*c) + 4*b*d**2*x**2*sqrt(-c**2*x**2 
+ 1)/(75*c**3) + 12*b*d*e*x**4*sqrt(-c**2*x**2 + 1)/(245*c**3) + 8*b*e**2* 
x**6*sqrt(-c**2*x**2 + 1)/(567*c**3) + 8*b*d**2*sqrt(-c**2*x**2 + 1)/(75*c 
**5) + 16*b*d*e*x**2*sqrt(-c**2*x**2 + 1)/(245*c**5) + 16*b*e**2*x**4*sqrt 
(-c**2*x**2 + 1)/(945*c**5) + 32*b*d*e*sqrt(-c**2*x**2 + 1)/(245*c**7) + 6 
4*b*e**2*x**2*sqrt(-c**2*x**2 + 1)/(2835*c**7) + 128*b*e**2*sqrt(-c**2*x** 
2 + 1)/(2835*c**9), Ne(c, 0)), (a*(d**2*x**5/5 + 2*d*e*x**7/7 + e**2*x**9/ 
9), True))
 
3.7.5.7 Maxima [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 314, normalized size of antiderivative = 1.30 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\frac {1}{9} \, a e^{2} x^{9} + \frac {2}{7} \, a d e x^{7} + \frac {1}{5} \, a d^{2} x^{5} + \frac {1}{75} \, {\left (15 \, x^{5} \arcsin \left (c x\right ) + {\left (\frac {3 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{2}} + \frac {4 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{4}} + \frac {8 \, \sqrt {-c^{2} x^{2} + 1}}{c^{6}}\right )} c\right )} b d^{2} + \frac {2}{245} \, {\left (35 \, x^{7} \arcsin \left (c x\right ) + {\left (\frac {5 \, \sqrt {-c^{2} x^{2} + 1} x^{6}}{c^{2}} + \frac {6 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{4}} + \frac {8 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{6}} + \frac {16 \, \sqrt {-c^{2} x^{2} + 1}}{c^{8}}\right )} c\right )} b d e + \frac {1}{2835} \, {\left (315 \, x^{9} \arcsin \left (c x\right ) + {\left (\frac {35 \, \sqrt {-c^{2} x^{2} + 1} x^{8}}{c^{2}} + \frac {40 \, \sqrt {-c^{2} x^{2} + 1} x^{6}}{c^{4}} + \frac {48 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{6}} + \frac {64 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{8}} + \frac {128 \, \sqrt {-c^{2} x^{2} + 1}}{c^{10}}\right )} c\right )} b e^{2} \]

input
integrate(x^4*(e*x^2+d)^2*(a+b*arcsin(c*x)),x, algorithm="maxima")
 
output
1/9*a*e^2*x^9 + 2/7*a*d*e*x^7 + 1/5*a*d^2*x^5 + 1/75*(15*x^5*arcsin(c*x) + 
 (3*sqrt(-c^2*x^2 + 1)*x^4/c^2 + 4*sqrt(-c^2*x^2 + 1)*x^2/c^4 + 8*sqrt(-c^ 
2*x^2 + 1)/c^6)*c)*b*d^2 + 2/245*(35*x^7*arcsin(c*x) + (5*sqrt(-c^2*x^2 + 
1)*x^6/c^2 + 6*sqrt(-c^2*x^2 + 1)*x^4/c^4 + 8*sqrt(-c^2*x^2 + 1)*x^2/c^6 + 
 16*sqrt(-c^2*x^2 + 1)/c^8)*c)*b*d*e + 1/2835*(315*x^9*arcsin(c*x) + (35*s 
qrt(-c^2*x^2 + 1)*x^8/c^2 + 40*sqrt(-c^2*x^2 + 1)*x^6/c^4 + 48*sqrt(-c^2*x 
^2 + 1)*x^4/c^6 + 64*sqrt(-c^2*x^2 + 1)*x^2/c^8 + 128*sqrt(-c^2*x^2 + 1)/c 
^10)*c)*b*e^2
 
3.7.5.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 598 vs. \(2 (215) = 430\).

Time = 0.31 (sec) , antiderivative size = 598, normalized size of antiderivative = 2.48 \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\frac {1}{9} \, a e^{2} x^{9} + \frac {2}{7} \, a d e x^{7} + \frac {1}{5} \, a d^{2} x^{5} + \frac {{\left (c^{2} x^{2} - 1\right )}^{2} b d^{2} x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )} b d^{2} x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{3} b d e x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {b d^{2} x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {6 \, {\left (c^{2} x^{2} - 1\right )}^{2} b d e x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {{\left (c^{2} x^{2} - 1\right )}^{4} b e^{2} x \arcsin \left (c x\right )}{9 \, c^{8}} + \frac {{\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b d^{2}}{25 \, c^{5}} + \frac {6 \, {\left (c^{2} x^{2} - 1\right )} b d e x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )}^{3} b e^{2} x \arcsin \left (c x\right )}{9 \, c^{8}} - \frac {2 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b d^{2}}{15 \, c^{5}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{3} \sqrt {-c^{2} x^{2} + 1} b d e}{49 \, c^{7}} + \frac {2 \, b d e x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{2} b e^{2} x \arcsin \left (c x\right )}{3 \, c^{8}} + \frac {\sqrt {-c^{2} x^{2} + 1} b d^{2}}{5 \, c^{5}} + \frac {6 \, {\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b d e}{35 \, c^{7}} + \frac {{\left (c^{2} x^{2} - 1\right )}^{4} \sqrt {-c^{2} x^{2} + 1} b e^{2}}{81 \, c^{9}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )} b e^{2} x \arcsin \left (c x\right )}{9 \, c^{8}} - \frac {2 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b d e}{7 \, c^{7}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )}^{3} \sqrt {-c^{2} x^{2} + 1} b e^{2}}{63 \, c^{9}} + \frac {b e^{2} x \arcsin \left (c x\right )}{9 \, c^{8}} + \frac {2 \, \sqrt {-c^{2} x^{2} + 1} b d e}{7 \, c^{7}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b e^{2}}{15 \, c^{9}} - \frac {4 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b e^{2}}{27 \, c^{9}} + \frac {\sqrt {-c^{2} x^{2} + 1} b e^{2}}{9 \, c^{9}} \]

input
integrate(x^4*(e*x^2+d)^2*(a+b*arcsin(c*x)),x, algorithm="giac")
 
output
1/9*a*e^2*x^9 + 2/7*a*d*e*x^7 + 1/5*a*d^2*x^5 + 1/5*(c^2*x^2 - 1)^2*b*d^2* 
x*arcsin(c*x)/c^4 + 2/5*(c^2*x^2 - 1)*b*d^2*x*arcsin(c*x)/c^4 + 2/7*(c^2*x 
^2 - 1)^3*b*d*e*x*arcsin(c*x)/c^6 + 1/5*b*d^2*x*arcsin(c*x)/c^4 + 6/7*(c^2 
*x^2 - 1)^2*b*d*e*x*arcsin(c*x)/c^6 + 1/9*(c^2*x^2 - 1)^4*b*e^2*x*arcsin(c 
*x)/c^8 + 1/25*(c^2*x^2 - 1)^2*sqrt(-c^2*x^2 + 1)*b*d^2/c^5 + 6/7*(c^2*x^2 
 - 1)*b*d*e*x*arcsin(c*x)/c^6 + 4/9*(c^2*x^2 - 1)^3*b*e^2*x*arcsin(c*x)/c^ 
8 - 2/15*(-c^2*x^2 + 1)^(3/2)*b*d^2/c^5 + 2/49*(c^2*x^2 - 1)^3*sqrt(-c^2*x 
^2 + 1)*b*d*e/c^7 + 2/7*b*d*e*x*arcsin(c*x)/c^6 + 2/3*(c^2*x^2 - 1)^2*b*e^ 
2*x*arcsin(c*x)/c^8 + 1/5*sqrt(-c^2*x^2 + 1)*b*d^2/c^5 + 6/35*(c^2*x^2 - 1 
)^2*sqrt(-c^2*x^2 + 1)*b*d*e/c^7 + 1/81*(c^2*x^2 - 1)^4*sqrt(-c^2*x^2 + 1) 
*b*e^2/c^9 + 4/9*(c^2*x^2 - 1)*b*e^2*x*arcsin(c*x)/c^8 - 2/7*(-c^2*x^2 + 1 
)^(3/2)*b*d*e/c^7 + 4/63*(c^2*x^2 - 1)^3*sqrt(-c^2*x^2 + 1)*b*e^2/c^9 + 1/ 
9*b*e^2*x*arcsin(c*x)/c^8 + 2/7*sqrt(-c^2*x^2 + 1)*b*d*e/c^7 + 2/15*(c^2*x 
^2 - 1)^2*sqrt(-c^2*x^2 + 1)*b*e^2/c^9 - 4/27*(-c^2*x^2 + 1)^(3/2)*b*e^2/c 
^9 + 1/9*sqrt(-c^2*x^2 + 1)*b*e^2/c^9
 
3.7.5.9 Mupad [F(-1)]

Timed out. \[ \int x^4 \left (d+e x^2\right )^2 (a+b \arcsin (c x)) \, dx=\int x^4\,\left (a+b\,\mathrm {asin}\left (c\,x\right )\right )\,{\left (e\,x^2+d\right )}^2 \,d x \]

input
int(x^4*(a + b*asin(c*x))*(d + e*x^2)^2,x)
 
output
int(x^4*(a + b*asin(c*x))*(d + e*x^2)^2, x)